Nonlocal interactions versus viscosity in turbulence 
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It is shown that nonlocal interactions determine energy spectrum in isotropic turbulence at small 
Reynolds numbers. It is also shown that for moderate Reynolds numbers the bottleneck effect is 
determined by the same nonlocal interactions. Role of the large and small scales covariance at the 
nonlocal interactions and in energy balance has been investigated. A possible hydrodynamic mech- 
anism of the nonlocal solution instability at large scales has been briefly discussed. A quantitative 
relationship between effective strain of the nonlocal interactions and viscosity has been found. All 
results are supported by comparison with the data of experiments and numerical simulations. 

PACS numbers: 47.27.-i, 47.27.Gs 



I. INTRODUCTION 



Classic wind-tunnel experiments [l[ showed that there is no scaling behavior in a presumably isotropic turbulence 
at low Reynolds numbers. Moreover, recent high-resolution numerical simulations [2j performed for low-Reynolds- 
number (R\ » 10 — 60) show no hint of scaling-like behavior of the velocity increments even when ESS [|| is applied. 
For moderate Reynolds numbers numerical simulations show 'excess' power just before the dissipation (Kolmogorov's) 
wavenumber kd (a hump in the compensated energy spectra), see for instance [J]. This non-scaling effect (bottleneck 
Q . effect @) is usually related with reducing efficiency of the energy cascade toward kd For the small Reynolds 

numbers applicability of the energy cascade idea is problematic for entire range of scales. It is shown in recent paper 
Q that nonlocal interactions become dominating in comparison with the local ones just in the near-dissipation range 
of scales (cf. also @). In this range the viscous effects cannot be neglected and scaling asymptote corresponding to 
I^h the nonlocal regime cannot be observed @ (in the inertial range the local interactions are presumably dominating 
ones and the nonlocal scaling asymptote also cannot be observed). However, the nonlocal scaling asymptote can be 
used as a zeroth term in a perturbation approach taking into account the viscosity effects. There are many ways 
. to develop such perturbation approach. For instance, in the paper Q a perturbation approach giving logarithmic 
corrections to the scaling was developed, provided by significant role of the kink instabilities of the vortex filaments at 
moderate and large values of Reynolds number fiol ]. This 'logarithmic'-perturbation approach is shown to be effective 
in a vicinity of the crossover scale r c , where exchange of stability between local and nonlocal regimes takes place 
Q. This vicinity is rather wide at moderate values of Reynolds number, when an overlap between these regimes is 
a strong phenomenon Q. For small Reynolds numbers, however, the kink instabilities of the small vortex tubes is 
suppressed by the strong viscosity (see, for instance [ll[ and references therein). Therefore, for the small Reynolds 
numbers another, adequate just for this case, perturbation approach should be developed. An approach of such type 
is suggested in present paper. Starting from this approach and using comparison with results of numerical simulations 
[l2| the laboratory experiments [l[ it is shown that energy spectrum for small Reynolds numbers is determined 
by the nonlocal interactions even in isotropic turbulence. The same nonlocal interactions provide a hydrodynamic 
mechanism for the so-called bottleneck effect for moderate Reynolds numbers. It is also shown that large and small 
scales covariance at the nonlocal interactions plays a significant role in these phenomena. A quantitative relationship 
between effective strain of the nonlocal interactions and viscosity has been found using dynamical equations. 



II. PERTURBATIONS TO SCALING 

Let us following to the paper || consider a dimensional function E(k) of a dimensional argument k. And let us 
construct a dimensionless function of the same argument 

If for kd 3> k we have no relevant fixed scale (scaling situation), then for these values of k the function a(k) must 
be independent on k, i.e. a(k) ~ const for kd S> k. For turbulence k could be a wavenumber and kd could be a 
dissipation wavenumber (kd — 1/?/, where n = (^ 3 /(e)) 1 / 4 is so-called viscous scale [24|). Solution of equation (1) 
with constant a can be readily found as 



E(k) ck 



(2) 
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FIG. 1: One-dimensional spectra measured in nearly isotropic turbulence downstream of a grid at small Reynolds numbers 
(the data are reported in [J]). The solid curves are drawn in the figure to indicate correspondence of the data to the equation 
(10) (non-local regime). 



where c is a dimension constant. This is the well-known power law corresponding to the scaling situations. 

Let us now consider an analytic approach, which allows us to find corrections of all orders to the approximate power 
law, related to the fixed scale kd- In the non-scaling situation let us denote 

/ = HE/A), x = \n(k/k d ) (3) 

where A and kd are dimensional constants used for normalization. 
In these variables, equation (1) can be rewritten as 

£ = a(*) (4). 

In the non-scaling situation a; is a dimensionless variable, hence the dimensionless function a(x) can be non-constant. 
Since the 'pure' scaling corresponds to k/k d C 1 we will use an analytic expansion in power series 

a(x) =a Q + ai(k/k d ) + ... + —a n {k/k d ) n + ... (5) 

where a„ are dimensionless constants. Choice of the small parameter for analytic perturbation approach is determined 
by physical situation, which one intends to consider. For instance, for moderate Reynolds numbers it would be generally 
preferable to consider the x~ x as a small parameter Q. In this paper, however, we intend to consider energy spectra 
in isotropic turbulence at small Reynolds numbers and corresponding phenomena in a relatively close vicinity of the 
dissipation scale at moderate Reynolds numbers. The kink instabilities of the vortex filaments, which are a significant 
factor near r c (see above and Q) are presumably not significant in situations with strong viscous effects [Tlj |. Therefore 
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FIG. 2: Three-dimensional spectra from a DNS performed in [l2l ] for different Reynolds numbers up to R\ = 125. The solid 
curves in this figure corresponds to the best fit by equation (10) (nonlocal interactions). The value of ai = 6.0 ± 0.1 is the 
same in all cases. 



the choice of x~ x does not seem to be relevant here. On the other hand, parameter k/kd seems to be less related 
to the specific hydrodynamic structures dominating processes in turbulence and more relevant to a sheer taking into 
account corrections to the scaling providing by viscosity. 

After substitution of the analytic expansion (5) into Eq. (4) the zcroth order approximation gives the power law (2) 
with a = OiQ. First order analytic approximation, when one takes only the two first terms in the analytic expansion 
(5), gives 

E(k) = ck ao e aik/kd . (6) 
(cf , for instance, [H|]-[l6j]). Corrections of the higher orders can be readily found in this perturbation approach. 



III. NONLOCAL INTERACTIONS 



Let us recall that in isotropic turbulence a complete separation of local and non-local interactions is possible in 
principle. It was shown by Kadomtsev [TtJ that this separation plays a crucial role for the local Kolmogorov's cascade 
regime with scaling energy spectrum 

E(k) ~ K (s) 2/3 k~ 5 / 3 (7) 

where (e) is the average of the energy dissipation rate, e, A; = 1/r is the wave-number, and K is the so-called 
Kolmogorov constant. This separation should be effective for the both ends. That is, if there exists a solution with 
the local scaling (7) as an asymptote, then there should also exist a solution with the non-local scaling asymptote. 
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FIG. 3: As in Fig. 2 but for R\ — 460. The upper part of the figure shows the same energy spectrum as the lower one but in 
the compensate (according to the Kolmogorov's scaling) form. The solid curves are the best fit corresponding to the nonlocal 
approximation (10) with the same «i = 6.0 ± 0.1 as in Figs 2. 



Of course, the two solutions with these asymptotes should be alternatively stable (unstable) in different regions of 
scales. It is expected, that the local (Kolmogorov's) solution is stable (i.e. statistically dominating) in inertial range 
(that means instability of the non-local solution in this range of scales). 

Roughly speaking, in non-local solution for small scales r only non-local interactions with large scales L(l>r/L) 
are dynamically significant (the interaction among the small scales is negligible compared with interaction via large 
scales) and the non-local interactions is determined by large scale strain/shear. This means that one should add to 
the energy flux (e)-parameter (which is a governing parameter for the both solutions) an additional parameter such 
as the strain s for the non-local solution. As far as we know it was noted for the first time by Nazarenko and Laval 
[HI that dimensional considerations applied to the non-local asymptote result in the power-law energy spectrum 

E(k) ~ c k- 1 (8) 
s 

both for two- and three-dimensional cases. Linear dependence of the spectrum (8) on (e) is determined by the linear 
nature of equations corresponding to the non-local asymptote that together with the dimensional considerations 
results in (8) (l8j . Interesting numerical simulations were performed in . In these simulations local and non-local 
interactions have been alternatively removed. For the first case a tendency toward a spectrum flatter than '-5/3' is 
observed near and beyond the separating scale (beyond which local interactions are ignored), that supports Eq. (8). 

Following to the perturbation approach suggested above both local and non-local regimes can be corrected. The 
first order correction is 

E(k) ~ K(e) 2 / 3 k- 5 / 3 e-^ k / k ^ (9) 

and 

E(k) si cMjfe- 1 e -ai(fc/fc ' t) (10) 
s 

for the local and non-local regimes respectively (K, c, a,\ and (3 are dimensionless constants). 
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Figure 1 shows one-dimensional spectra measured in nearly isotropic turbulence downstream of a grid at small 
Reynolds numbers (the data are reported in The solid curves are drawn in the figure to indicate correspondence 
of the data to the equation (10) (non-local regime). The experimental data, however, is not controlled enough to 
study fine properties of the isotropic turbulence. Therefore below we will mainly use the data obtained in numerical 
simulations. 

It is shown in that there is an 'exchange of stability' phenomenon at certain k c = l/r c . That is, for k < k c the 
Kolmogorov's regime is stable and the non-local regime is unstable, whereas for k > k c the Kolmogorov's regime is 
unstable and the non-local regime is stable. For this scenario, at k — k c the Kolmogorov's regime is still asymptotically 
scale-invariant (i.e. Eq. (7) gives an adequate approximation for this regime), while for the non-local regime the first 
order correction is substantial (i.e. Eq. (10) should be used at k > k c for the non-local regime). In this scenario the 
Kolmogorov's regime plays significant role in the viscous stabilization of the non-local regime for k > k c (see [IH), but 
for these k the non-local regime becomes statistically dominating instead of the Kolmogorov's one (it is well known 
that only for energy spectra steeper than fc~ 3 it can be rigorously proved that the dominant interactions are nonlocal 
(cf Eq. (10) that provides such steepness). 

It should be noted that in [20j a second scaling solution E(k) oc Pk^ 1 was obtained in addition to the Kolmogorov 
scaling using the Clebsch formulation of hydrodynamics. The prefactor P — ((e)i/) 1 / 2 denotes a certain flux in the 
wavenumber space [20]. It will be shown below that the prefactor in Eqs. (8), (10) can be transformed into the one 
formally equal to P (see Eq. (15) and Fig. 4). This could mean that the studied nonlocal regime (more precisely: 
its scaling asymptote) is closely related to the Yakhot-Zakharov scaling solution introduced in the Ref. [20j. Such 
identification (if valid) allows us understand instability of the nonlocal regime for large scales (small k). Indeed, it is 
shown in [2l[ that reconnection process breaks conservation of the integral determining the flux P for the large scales 
(small k) and , therefore, there is no possibility for realization of the scaling asymptote k~ l of this solution. Following 
to Newell [2l| [22j the nonconservation of the integral in large scales generally follows from nonlocality of the viscosity 
term in the equations formulated for the Clebsch variables. Since for the large scales the nonlocal regime should be 
represented by its asymptote E oc k^ 1 (as an intermediate asymptote [23| the nonconservation of the integral makes 
this regime unstable for the large scales (small k). 

Figure 2 shows three-dimensional energy spectra calculated using data from a high-resolution direct numerical 
simulation of homogeneous steady three-dimensional turbulence [ijj for different Reynold numbers up to R\ = 125. 
The solid curves in this figure corresponds to the best fit by equation (10) (nonlocal interactions). The exponent 
a\ ~ 6.0±0.1 for all considered values oiR\. The arrows indicate the r c = l/k c scales calculated using corresponding 
D LLL (r) (see @). 

In figure 3 we show the data obtained in the same DNS as those shown in Fig. 2 but for R\ — 460. In the upper 
part of this figure we show the energy spectrum in the compensate (according to Kolmogorov's scaling Eq. (7)) form. 
One can clear see the hump corresponding to the bottleneck effect. The solid curves are the best fit to the the same 
nonlocal spectrum (10) with the same (universal) value of a\ ~ 6.0 ± 0.1 (see also below). One can see that the 
bottleneck effect is determined by just the same nonlocal interactions as the above considered energy spectra for small 
Reynolds numbers. 



IV. STRAIN AND VISCOSITY 

It is shown in Ref. Q that for sufficiently large Reynolds numbers, providing a visible inertial interval, there is 
an overlapping between the two regimes: non-local and local (Kolmogorov). This overlapping is based on the very 
nature of the stability exchange between the two statistical regimes. However, in analogy with the viscous scale 
7] = (v 3 / (s)) 1 / 4 [11] there should be a scale rj s such that for scales r < r) s contribution of the local interactions will 
be drastically decreased in comparison with the nonlocal ones. In the analogy with rj one can calculate rj s using the 
dimensional considerations as 

(cf for shear flows [25|). 

The dynamical equations provide us with a relationship [24] 



(e) = 1v I k 2 E{k)dk (12) 
Jo 

For E(k) given by Eq. (10) the dissipation function k 2 E(k) has its maximum at k = fc^/6. If kd 3> k s — l/rj s (see 
below) the maximum of the dissipation function k 2 E(k) is located just between k s and kd- Therefore, we can estimate 
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(<s>v) 1/2 

FIG. 4: The prefactor in the nonlocal approximation to the energy spectra (10): c(e)/s - circles, versus {{e)u) 1/2 . The straight 
line with the slope equals to 18 indicates agreement with Eq. (15). 



(12) as 

(e) ~ 2v / k 2 E{k)dk (13) 

Jk s 

Substituting (10) (with a% — 6) into (13) we obtain relationship 



2 , 



s^± { e)^u-^ (14) 
Using the relationship (14) one can estimate the prefactor in the nonlocal approximation to the energy spectra (10) 

as 

cM ~ 18((eM 1/2 (15) 
s 

Now using the data of the DNS [l2| (cf Figs. 2,3) let us calculate the prefactor c(e)/s. Results of these calculations 
are shown in figure 4 as circles (Rc\ — 38,54,70,125, 284,380,460). This figure shows the prefactor c (e)/s against 
((e)^) 1 / 2 . The straight line with the slope equals to 18 indicates agreement with Eq. (15). 

It should be noted that in the DNS [12| (e) — const for R\ > 70, in agreement with the well known Kolmogov's 
hypothesis [24|. Therefore, for R\ > 70 (when (e) ~ const [12]) it follows from Eq. (14) that s oc v^ 1 / 2 . This 
relationship provides us also with dependence of the strain s on R\. 

Due to the nonlinear character of the Navier-Stokes equations the so-called triadic type of interactions is dominating 
mechanism of the dynamical interactions in turbulence (see, for instance, Q). A triad corresponding to the nonlocal 
interactions involves two short- wave- number modes and one long- wave-number mode (see a sketch in figure 5). Ac- 
cordingly, at the nonlocal interactions two characteristic space scales are actively involved: large-scale characteristic 
scale rj s (11) and small-scale (viscous or Kolmogorov) characteristic space scale n = (^ 3 /(e)) 1 / 4 . It is naturally that 
the large scales should be normalized by n s while the small scales should be normalized by n. This could cause 
an obvious problem at the nonlocal interactions. However, the large and small scales covariance relationship at the 
nonlocal interactions follows directly from the relationship (14) 



/ / \ \ 3/4 / 1 o \ 3/2 



const (16) 



(i.e. the relation rj s /rj is independent on R\). And vice versa, the large and small scales covariance (15) results in the 
relationship of the type (14). The covariance at the nonlocal interactions supports right balance between the energy 
flux to the small scales and their dissipative capacity (cf, for instance, [26| . [27| ) . 
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FIG. 5: Sketch of the nonlocal triadic interaction in the physical and in the Fourier space 

Actually, the scale Ws = l/k s should provide an edge of applicability of the approximation (10) to the real energy 
spectra. If one compares the rj s determined by this way from the Figs. 2,3,6 one can see that indeed ?7 s /?7 — const 
in agreement with Eq. (16) (we have also checked this with the data for R\ = 284,380 [12] and the data reported 
in [29(). Moreover, taking into account a continuity condition of the energy spectrum at the point k = k s (and using 
Eqs. (7) and (10)): 

K (e) 2 / 3 k^ 3 ~ c M fc -i e -«i(WM n 7) 
s 

we obtain equation 

K ~ c exp[-(2c) 3/2 /a?] (18) 

Substituting the Kolmogorov constant K ~ 1.6 (see [Hj],[28j]) and ot\ ~ 6 into this equation we obtain c ~ 2. Then, 
substituting this value of c into Eq. (16) we obtain k s /kd — 0.037. The last value is in agreement with the available 
data (see Figs. 2,3,6 and (H). 

Together with the universality of cki — 6.0 (which also is a consequence of the scale covariance) Eq. (16) determines 
the well known from DNSs (l2l|.(29l| universality (independence on R\) of the position of the 'hump' in the axes kr\ 
for the bottleneck effect. 



V. PASSIVE SCALAR 

For passive scalar 9 in the isotropic turbulence the equations (9), (10) should be replaced Q by the equations 

E 9 {k) oc {e)- 1/3 {e e ) k~ 5/3 e-^ k / k ^ (19) 

and 

E \k) oc M fc -i e -Wfe«) (20) 
s 

for the local and non-local regimes respectively (7 and S are dimensionless constants, and (eg) is the average value of 
dissipation rate of scalar variance). 

Figure 6 shows three-dimensional passive scalar spectrum from a DNS performed in [30j for Peclet number P\ = 427 
(the Schmidt number is unity, i.e P\ — R\). The solid curves in this figure correspond to the best fit by equation (20) 
(nonlocal interactions). In the upper part of this figure we show the spectrum in the compensate (according to the 
Corrsin-Obukhov scaling [3(|) form. One can clear see the hump corresponding to the bottleneck effect. The arrows 
show position of the rj s scale. Comparing with Figs. 2,3 one can see that for the passive scalar r) s /r] takes the same 
universal value as for the velocity field (the large and small scales covariance (15)). 
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FIG. 6: The passive scalar DNS data (circles) for homogeneous isotropic turbulence described in [30|], Peclet number P\ = 427 
and the Schmidt number is unity (i.e Pa = R\)- The solid curves in this figure correspond to the best fit by equation (20) 
(nonlocal interactions) . 
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